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ABSTRACT 


A  semi -empirical  analytic  model  for  leeside  flow  field  behavior  over 
sharp  circular  cones  at  large  incidence  is  described.  The  lJietbodology 
volves  replacement  of  the  actual  body  by  an  mvlscld  equivalent  body  over 
which  the  flow  field  Is  developed.  The  equivalent  body  Is  selected  by  a 
combination  of  turbulent  boundary  layer  analysis  up  to  the  primary  separa 
tion  point  and  correlation  of  available  experimental  results  within  the 
separated  regions.  The  inviscid  flow  field  over  the  equivalent  body  I s 
developed  by  means  of  a  numerical  finite  difference  marching  technique 

(Kutler  program)  . 

Computer  programs  for  implementing  the  analysis  were  developed  and 
delivered  to  the  AFFDL/FXG  together  with  user's  manuals.  The  performance 
of  these  programs  is  demonstrated  by  several  numerical  examples  and  com¬ 
parisons  with  experiment. 
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SECT  I  ON  I 


INTRODUCTION  AMD  BACKGROUND 


Up  until  the  last  decade  very  little  interest  o  *-  te  1  •"  deve  i  -m- 
of  computational  methods  for  describing  flow  *  '  e  1  c  ben, sc  I  _>r  c  Ices 

of  conical  bodies  at  very  large  incidence,  This  i  ac*  of  .  ntcres  t  was. 
consequence  primarily  of  practical  considerations  since,  unt'l  that  1  i R  , 
sustained  flight  at  angles  of  attack  comparable  to  o-  greater  than  the  .w 
hide  thickness  to  length  ratios  was  not  required  to  fulfill  any  pract  :: 
mission  requirements. 

With  the  advent  of  the  space  shuttle  corcept  and  tie  rapid  improv  crw- - 
of  ABM  performance  this  situation  has  changed  dramatically,  as  evidenced 
the  many  studies  of  this  problem  undertaxen  during  the  post  several  yu j 
As  might  be  expected,  many  of  these  studies  were  of  an  exper ' nente I  rat 
(References  1-6),  with  most  analysis  restricted  to  treatment  of  cones  as 
angles  only  up  to  about  the  cone  haif  angle  (References  7" 1 3).  The  problem 
of  course,  is  an  exceedingly  difficult  one  since  it  involves  the  growt' 
separation  of  a  three  dimensional  boundary  layer  wh  i  ch  subsequent  ’  *  >  n  U 
with  the  outer  inviscid  flow.  This  phenomenon  can  effect  vehicle  perfa 
ae  rodynam  i  ca  I  1  y  (stability)  as  well  as  thermodynami  ca  I  1  y  (heating  rate-,  1 
location  of  the  separation  lines  and  the  strength  of  the  vortices  wr 
from  rollup  of  the  shear  layer  evidently  can  influence  the  flight  cIm"  : 

a  given  vehicle  configuration.  Furthermore,  with  regard  to  the  heating 
tion,  it  was  found  »xper i mental !y  that  heat  transfer  rates  on  the  let^ar- 
faces  could  achieve  levels  comparable  to  those  on  the  windward  side  and,  r 
certain  conditions,  even  exceed  these  values.  The  practical  implication 
these  observations  was  apparent  and  provided  the  motivation  for  develop! 
analytic  techniques  capable  of  describing  the  details  of  the  phenomenon 

Prior  to  ATI's  current  study,  such  analyses  as  were  available  were 
stricted  to  reiatively  low  angles  of  attack  and  laminar  conditions.  Be*' 

(13)  is  a  representative  example  of  these  studies  but  the  flow  fieln  .rr;  e-1 
therein  differs  in  a  fundamental  way  from  that  investigated  by  ATL  dunir  tie 
past  year. 
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In  order  to  clarify  these  differences  it  is  instructive  to  examine  and 
compare  the  information  presented  in  Figures  (la),  (lb)  and  (lc)  with  tha*  in 
Figures  (2a),  (2b)  and  (2c).  It  is  evident  from  this  comparison  that  there 
are  two  distinct  flow  patterns  associated  with  the  small  and  large  angles  of 
attack.  The  first  type  consists  of  one  separation  point  and  a  single  re¬ 
circulating  region.  The  other  consists  of  two  separation  points  and  two  re¬ 
circulation  regions.  It  has  been  concluded  from  our  research  effort  during 
the  past  year  that  the  single  separation  exists  for  the  following  conditions: 

1.  Low  angles  of  attack  (a/9^  <  1.2)  for  both  laminar  and  turbulent  con¬ 
ditions,  high  and  low  Mach  numbers. 

2.  High  angles  of  attack,  M^  >  5,  and  low  Reynolds  number. 

The  double  separation  type  of  flow  field  which  has  been  investigated  at  ATL 
exists  for  the  following  conditions: 

(a)  High  Reynolds  number  (turbulent)  for  all  ranges  of  Maci  number 
and  for  ot/6^  >  1.2. 

(b)  Low  Mach  number,  laminar  flow  conditions  and  a/6c  >  1.2. 

A  direct  comparison  of  these  two  patterns  is  depicted  schematically  in  Figure 
(3).  The  associated  pressure  distributions  are  shown  in  Figure  (L)  .  While 
conditions  for  laminar  flow  could  be  analytically  treated  and  hfjs  been  treated 
recently  for  the  single  type  of  separation,  the  double  separation  flow  field 
for  turbulent  conditions  is  complex  and  requires  a  semi -emp i r i cal  type  of 
anal ys i s . 

From  a  practical  point  of  view,  the  motivation  for  developing  an  analytic 
capability  for  describing  the  second  type  of  flow  phenomenon  ultimately  derives 
from  a  need  to  compute  heat  transfer  characteristics  on  the  leeward  surfaces 
when  these  double  recirculation  zones  or  vortices  are  present.  Note  that  con¬ 
siderable  evidence  exists  which  indicates  that  the  observed  high  heating  rates 
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FIGURE  2b  TYPICAL  PRESSURE  DISTRIBUTION  FOR  THE  DOUBLE 
SEPARATION  CASE-COMPARISON  OF  THEORY  AND 
EXPERIMENT 
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FIGURE  2C  CONCEPTUAL  REPRESENTATION  OF  FLOW  FIELD  FOR  THE  DOUBLE 
SEPARATION  CASE 
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mentioned  earlier  are  directly  attributable  to  the  formation  or  such  vortices. 


Under  certain  conditions  associated  with  practical  flight  trajectories  of 
blunted  cones,  regions  of  laminar  and  turbulent  conditions  can  occur  simulta- 

Ineously.  It  is  conceivable  that  this  transition  can  trigger  the  vortices  which 
would  not  occur  under  purely  laminar  conditions.  Evidently,  the  availability 
of  a  reliable  flow  field  predictive  technique  is  a  prerequisite  for  determina¬ 
tion  of  when  and  if  these  mechanisms  occur. 


In  some  recent  experiments  conducted  at  the  New  York  University  Aerospace 
Laboratory  (Reference  14),  some  peculiarities  with  regard  to  heat  transfer 
measurements  have  been  observed.  These  could  be  explained  in  the  light  of 
the  phenomena  described  above.  A  summary  of  these  results  is  presented  in 
Figure  (5).  In  this  case,  heat  transfer  mf-'-.  on  the  leeward  plane  are  not 
only  higher  than  those  on  the  windward  aminar  conditions  but 

under  turbulent  conditions  as  well.  none  associated  with  these 

regions  of  high  heating  may  be  seer  cture  (Figure  6)  where 

the  feather  type  markings  are  locat< 


Similar  behavior  has  been  observed  shuttle  type  vehicles  by 

Whitehead  at  Langley  Field  (Reference  15).  Figure  (7)  shows  the  dependence 
of  these  heating  rates  on  Reynolds  number  while  Figure  (8)  shows  the  variation 
with  Reynolds  number,  Mach  number  and  angle  of  attack.  It  is  evident  from 
these  results  that  vortex  heating  is  significant  and  needs  to  be  accounted  for 
to  insure  proper  design  of  high  performance  reentry  configurations. 

The  method  of  approach  and  analysis  utilized  by  ATL  to  treat  this  problem 
during  the  past  year  is  described  in  the  next  two  sections.  Results  and  com¬ 
parisons  with  experiment  are  presented  in  the  subsequent  section. 
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FIGURE  5  MEASURED  HEA1  TRANSFER  RATES  ON  THE  SURFACE  OF  BLUNTED  CONES  AT  INCIDENCE-DATA  FROM  REFERENC 
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FIGURE  7  EFFECT  OF  REYNOLDS  NUMBER  ON  CENTERLINE  LEE-SURFACE  HEATING 


FIGURE  8  HEAT  TRANSFER  COEFFICIENT  ON  LEE  SIDE  OF  CIRCULAR  CONES 


SECTION  I  I 

METHOD  OF  APPROACH  AND  ANALYSIS 


The  overall  objective  of  the  current  effort  was  to  develop  a  numerical 
computation  scheme  which  would  provide  a  description  of  the  flow  field  char¬ 
acteristics  and  surface  conditions  on  the  leeside  of  cones  at  large  incidence 
and  at  conditions  corresponding  to  a  turbulent  boundary  layer  state.  The 
method  of  approach  utilized  to  achieve  this  objective  is  a  semi-empirlcal  one 
which  utilizes  an  inviscld  equivalent  body  in  order  to  develop  the  flow  field. 
The  equival  nt  body  is  derived  from  a  combination  of  viscous  turbulent  bound¬ 
ary  layer  program  up  to  the  primary  separation  point  and  a  correlation  of  flow 
field  data  within  the  separated  region.  The  method  involved  the  following  in¬ 
gredients  : 


A*  Physical  Model  -  The  physical  model  which  was  employed  incorporates 
the  concept  of  an  "effective"  body  to  represent  the  region  of  high  shear  ad¬ 
jacent  to  the  geometric  body.  That  is,  the  shear  layer,  Including  both  the 
boundary  layer  proper  and  the  vortices,  was  cons-idered  to  be  enclosed  by  a 
streamline  or  stream  surface  which  divided  the  region  of  high  shear  from  the 
outer  invlscid  shock  layer.  The  pressure  across  this  shear  layer  is  assumed 
constant  and  the  invlscid  flow  field  develops  as  if  in  the  presence  of  a  body 
whose  cross  sectional  geometry  Is  defined  by  the  aforementioned  dividing  stream 
surface. 


B*  Flow  Field  Computation  Scheme  -  Once  the  effective  body  shape  has  been 
determined  a  computation  scheme  capable  of  determining  a  three  dimensional, 
supersonic  flow  field  Is  required.  Actually,  as  will  be  discussed  below,  this 
scheme  is  also  utilized  to  aid  in  the  determination  of  the  effective  body 
shape.  The  flow  field  computation  scheme  (FFCS)  selected  for  this  purpose  is 
that  due  to  Kutler,  Reference  (16)  which  was  specifically  designed  to  analyze 
the  flow  field  over  three  dimensional  configurations  at  high  angles  of  attack 
at  NASA-Ames. 

The  analysis  involves  numerical  integration  of  the  finite  difference 
analogue  of  the  conservation  equations  for  three  dimensional  invlscid  iso- 
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energetic  flow.  The  working  independent  variables  in  this  numerical  scheme 
are  related  to  the  cylindrical  coordinates  z,  r,  and  $  according  to 


5(z,r,*)  -  (r-rb)/(rs"rb) 

nU)  -  arctan  (k  tan  4>) 

where  z  is  aligned  with  the  body  axis.  Here  rb*=  r^ (z , <}>)  and  rs**r5  (2  represent 
the  vehicle  body  geometry  and  shock  location  respectively.  The  parameter  k  Is 
a  free  parameter  which  serves  the  purpose  of  distorting  the  mesh  in  the  physical 
plane;  i .  e .  :  equal  increments  of  n  results  In  a  clustering  of  cj>  -  i  ncrements  in 
the  vicinity  of  <j>-90°  when  0  <  k  <  1  while  for  k  >1  the  clustering  occurs  in 
the  vicinity  <(>-0,  180°.  For  <  >  1  of  course  the  spacing  approaches  uniformity. 


The  integration  scheme  involves  a  straight  marching  procedure  In  the 
hyperbolic  coordinate  z.  Accordingly,  the  scheme  can  be  applied  only  to  a 
supersonic  domain  and  further  requires  specification  of  valid  initial  data 
along  some  starting  line.  For  the  general  case  this  must  be  supplied  by  other 
means  (e.g.:  a  subsonic  blunt  body  scheme).  For  the  present  application,  how¬ 
ever,  where  body  geometries  are  restricted  to  the  conical  case  (dr^/dz  ■>  constant) 
initial  data  can  be  arbitrary  since,  for  sufficiently  large  z,  the  desired 
conical  solution  will  be  obtained  independent  of  the  correctness  of  the  initial 
data. 


Details  of  the  analysis  and  the  numerical  methods  associated  with  the  com¬ 
putation  scheme  are  not  included  here  since  they  are  well  documented  in  Refer¬ 
ence  (16).  In  particular,  it  is  noted  that  ATL's  involvement  in  the  develop¬ 
ment  of  the  FFCS  consisted  solely  of  developing  a  specialized  geometry  package 
suitable  for  the  present  purpose.  This  geometry  package  has  been  substituted 
for  the  version  incorporated  by  the  original  computer  code  as  received  from 
NASA-Ames. 
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The  necessity  for  revising  the  geometric  representation  is  exemplified  by 
the  comparison  shown  in  Figure  (9).  As  can  be  seen,  the  general  character¬ 
istics  of  the  geometry  needed  to  represent  the  shear  layers  on  cones  at  high 
incidence  differ  in  a  very  profound  way  from  the  shuttle-like  configuration 
for  which  the  original  computer  code  was  intended. 

The  geometry  used  to  represent  the  shuttle  configuration  in  Reference 
(16)  is  indicated  in  Figure  (10)  and  is  labelled  Kutler's  geometry  where  it 
is  compared  with  our  initial  proposed  modifications.  The  proposed  changes 
consisted  of  the  addition  of  two  sets  of  four  segments  between  straight  line 
segment  5~t>  and  th°.  original  ellipse  6-7*  Each  set  of  segments  consisted  of 
ellipse,  line,  cubic  and  line.  The  additional  geometrical  segments  were  as 
follows:  ellipse  -  segment  6-7  (defined  using  point  i),  line  -  segment  7-8, 
cubic  -  segment  8-9,  line  -  segment  9"10,  ellipse  -  segment  10-11  (defined 
using  point  j),  line  -  segment  11-12,  cubic  -  segment  12-13  and  line  -  segment 
13"!^.  The  original  ellipse  6-7  was  replaced  by  the  equivalent  ellipse  14-15. 
Specification  of  these  8  additional  segments  required  16  input  parameters 
namely  Xu,  Y^,  Y Yg,  X^,  Y^,  0^,  Y)Q,  Xh>  Y^,  Y^,  Y]2,  Xjg,  Y^,  0^  and 
Y^  with  Yj^  replacing  Y^  of  the  original  input.  Figure  (11)  is  a  chart 
indicating  the  changes  which  were  required  to  the  input  of  the  program. 

These  modifications  were  implemented  and  a  series  of  calculations  were 
attempted  with  this  version  of  the  computer  code.  It  was  found,  however,  that 
even  this  augmented  version  of  1A  segments  was  unsatisfactory  in  terms  of 
providing  the  flexibility  to  reproduce  the  desired  pressure  variations.  Ac¬ 
cordingly,  a  further  modification  was  implemented  involving  a  straightforward 
tabular  input  of  effective  body  radius  at  each  of  the  peripheral  grid  points. 
This  scheme  proved  to  be  successful  and  determination  of  effective  body  shapes 
was  achieved  without  further  difficulty. 

Once  this  had  been  accomplished  these  body  shapes  were  generalized  in 
terms  of  their  dependence  on  displacement  thickness  at  the  primary  separation 
point  and  angle  of  attack  as  is  discussed  in  Section  I  I  I E.  The  logic  of  the 
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FIGURE  10  SCHEMATIC  OF  KU1LER  AND  PROPOSED  GEOMETRY  FOR  BODIES  AT  ANGLE 
OF  ATTACK 
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Figure  11,  INPUT  FOR  GEOMETRY  DESCRIPTION 
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geometry  subroutine  was  subsequently  modified  so  that  the  requisite  tabular 
inputs  are  generated  internally  consistent  with  this  generalization.  This 
aspect  is  discussed  in  detail  in  Section  I  1  I E ,  as  well  as  in  the  user's  manual 
for  the  FFCS  which  has  been  delivered  to  the  AFFDL/FXG. 

C.  Boundary  Layer  Computation  Scheme  -  A  boundary  layer  computation 
scheme  (BLCS)  is  needed  ror  two  reasons.  First,  in  order  to  characterize 
the  effective  body  shape  in  the  separated  region  the  boundary  layer  state 
up  to  the  primary  separation  point  needs  to  be  defined.  Furthermore,  once 
the  effective  body  shape  has  been  determined,  boundary  layer  calculations 
within  this  region  are  needed  to  define  the  surface  conditions. 

The  boundary  layer  scheme  selected  for  this  purpose  utilizes  the  stream¬ 
line  tracing  concept  with  a  small  cross-flow  assumption  proposed  by  Vagi  f  o - 
Laurin  References  (17)  and  (18).  Details  of  the  analysis  may  be  'ound  in 
References  (19)  and  (20).  These  also  describe  the  numerical  procedures  and 
computer  code  developed  to  implement  this  scheme.  A  user's  manual  for  the 
BLCS  Is  also  available  from  the  AFFDL/FXG.  The  selected  method  provides 
the  three  dimensional  capability  which  is  an  essential  feature  of  the  viscous 
flow  phenomena  involved  in  the  present  study. 

It  should  be  noted  that  implementation  of  this  computation  scheme  requires 
a  priori  knowledge  of  the  inviscid  pressure*  variation  along  streamlines  as 
well  as  the  distribution  of  h^.  the  parameter  which  characterizes  the  stream¬ 
line  spreading.  Determination  of  these  properties  requires  a  streamline 
tracing  capability.  During  the  current  contractual  effort  a  special  scheme 
valid  for  sharp  nosed  circular  cones  was  developed  which  permitted  boundary 
layer  calculations  to  be  performed  up  to  the  primary  separation  point  only. 

A  description  of  this  analysis  is  given  in  the  next  section. 

*For  blunted  cones  the  variation  of  edce  entropy  is  also  needed.  This  applica¬ 
tion  has  not  been  considered  during  the  present  study. 
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D.  Streamline  Tracing  Analysis  -  The  assumptions  on  which  this  analysis 
is  based  include  that  the  flow  is  conical,  the  bow  shock  is  straight  and  at¬ 
tached  and  the  flow  downstream  of  the  shock  is  isentropic.  Under  these  cir¬ 
cumstances  the  components  of  velocity  on  the  cone  surface  in  the  peripheral 
direction  v^  is  related  to  that  along  cone  generators  v^  by  the  differential 
equation 


d  VR 

-r r “  “  V.  sin  0 
d<p  <(>  c 


(4) 


where  =  v^/V^  and  =  v^/V^  are  the  velocity  components  normalized  with 
respect  to  the  limiting  velocity  and  0  is  the  cone  half  angle.  The  coordi¬ 
nate  notation  is  defined  in  Figure  (12).  Equation  (1)  follows  from  Equation 
(14-4)  of  Reference  (21)  when  applied  to  the  surface  of  the  cone  . 


The  velocity  components  can  also  be  related  to  the  local  pressure  via  the 
equat i on 


£_ 

Pco 


00 


exp (-AS/R) 


where 

2  2  2 
V  +  < 


(5) 

(6) 


V(M  is  the  free  stream  velocity  normalized  with  respect  to  VL  and  AS/R  represents 
entropy  jump  across  the  oblique  shock  on  the  windward  side  in  the  plane  of 
e  t  ry ;  i .  e .  : 


exp(-AS/R)  =  p  /p  (7) 

2 

The  quantity  AS/R  is  a  constant  under  the  assumed  conditions  and  can  be  evaluated 
once  and  for  all  when  the  shock  incl incation  is  known. 
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If  the  pressure  distribution  is  prescribed  as  a  function  of  $  Equations 
W  *  (5)  and  (6)  represent  three  equations  for  the  determination  of  the  two 
velocity  components  and  VR  and  the  resultant  V.  This  is  the  approach 
adopted  here.  Specifically,  the  pressure  variation  is  considered  to  be  given 
either  as  a  tabular  input,  or,  alternately  is  represented  by  curve  fitting 
three  points  in  accordance  with  the  trigonometric  relation 

P(4>)  “A  cos  <j>  +  B  cos  2(j>  +  C  (8) 

where 

P  5  p/p 

The  coefficients  appearing  in  Equation  (8)  are  evaluated  by  specifying  the 
pressures  Pj,  P^  and  P^  at  values  of  the  peripheral  angles  <j>2  and  <(> 
respectively.  Then 

A  -  [P,(B2-B3)  +  P2(B3-B1)  +  P3(B,-B2)]  d’1  (9a) 

b  -  (p,(a3-a2)  +  p2(a,-a  )  +  p3(a2-a1)i  d'1  (9b) 

C  “  [P1(A2B3“A3B2)  +  P2(A3B1_A1B3)  +  p3(a,B2-A2B))]  d"1  (9c) 

where 

o  -  a2b3*a3b,  ♦  a,b2  -  a28,  -  a,b3  -  a3b2  Od) 

and 

Aj  =  cos  <j>j 
Bj  =  cos  2  4>  j 

Equation  (8)  has  been  found  to  provide  a  good  representation  of  pressure 
variation  around  sharp  circular  cones  at  incidence  (Reference  22)  at  least 
for  angles  0  K  $  <  120°. 
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With  P(<(>)  specified  the  variation  of  the  resultant  velocity  V($)  can  be 

obtained  in  an  elementary  manner  from  Equation  (5).  Determi  nat  ion  of  V  (<J>) 

$ 

and  VRU)  requires  numerical  integration  of  Equation  (4).  The  procedure 
utilized  here  is  as  follows.  In  terms  of  the  finite  difference  mesh  defined 
in  Figure  (13)  we  can  write  the  central  difference  form  of  Equation  (4)  as 

VRU)  ■  V^(i-2)  +  2  V  ( i - 1 )  A<j>  sin  @c  (10) 

Equation  (10)  permits  calculation  of  VD  at  a  generic  point  i  within  an  error 
2  * 

of  order  ( Acp )  in  terms  of  data  evaluated  at  the  two  preceding  mesh  points. 

The  corresponding  value  of  is  given  by  virtue  of  Equation  (6)  as 

%(0  -  [V2  ( i )  -  vjj(i)] 

where  we  recall  that  V(i)  is  known  at  all  <j>(i).  To  initiate  this  procedure 
we  apply  the  boundary  condition 

y>  ■  ° 

VR(I)  -  V(l) 

which  follows  from  symmetry  considerations  since  )  =  0.  We  also  make  use 
of  the  one  sided  difference  form  of  Equation  (4);  viz. 

V,(2) 

VR(2)  -  VR(l)  +  “2 —  4<j>  sin  6c 

also 

V2 (2)  +  V2( 2)  -  V2 (2) 

Elimination  of  VR(2)  between  these  two  equations  yields  the  quadratic  relation 

a  V2 (2)  +  b  V  (2)  +  c  =  0  (12) 


(I  la) 

(lib) 
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where 

A4>  sin  0  2 

a  -  1  +  ( - 5 - £■) 

b  *  VD(l)  *A<J>  sin  0 

K  C 

c  -  y2(i)  -  v2 (2) 

Equation  (lib)  and  (12)  permit  evaluation  of  V,(2)  and  V_ (2) ,  in  this  case, 

$  K 

within  an  error  of  order  A$.  Accordingly  V.,  V_  and  V  can  be  determined  for 

9  K 

all  <J>  by  a  step  by  step  integration  procedure  using  Equation  (10).  Certain 
auxiliary  relations  which  are  of  interest  can  then  also  be  computed.  In 
particular,  the  flow  inclination  relative  to  cone  generators  b  Is  determined 
from  the  relation 


ip  -  arc  tan 


(13) 


The  derivative  3V^/3b  is  also  needed.  Following  Reference  (22)  this  is 


given  by 


i  dA 

7  3<t> 


•  2 
sin  e 


where 


2V 


(y-i)  (\-r) 


i 


-JL  (i  lZ)i 

2  'P  2/J 
y  M  30 

e 


sin  0 


(14) 


(15) 


For  the  case  when  the  pressure  distribution  is  curve  fit,  the  second 
derivative  of  pressure  follows  immediately  from  Equation  (8);  viz: 


32P 

3$2 


(A  cos  <(>  +  4  B  cos  24> ) 


06) 


For  the  case  where  the  pressure  distribution  is  represented  by  a  tabular 
input  of  P  vs.  $  Equation  (14)  is  replaced  by 


(i  i!t) 

Vv  3b  ’ 


v  (i+D  -  v  (i ) 
■"•'vTiTIb  — 


(17) 
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With  the  conical  flow  field  thus  established  tracing  of  streamlines  on 
the  surface  of  the  cone  requires  determination  of  the  <=  treamwise  increment  ds 
associated  with  an  arbitrary  increment  dR  along  a  conical  ray  starting  at  any 
generic  point"  Q(Rj ,  *  ).  This  is  given  by 

2  2  2  2 

ds  -  dR  +  r l  d,2  (I8) 

where  rc  is  the  local  radius  of  the  cone  which  is  a  known  function  of  R;  t.e.: 
rc  "  R  s'n  ®c*  i n  Equation  (18)  db  is  obtained  from 

tan  \p.  dR 

d4>  "  R.  +JdR  (19) 

which  follows  by  definition  from  Equation  (13).  |n  Equation  (19)  the  stream¬ 
line  inclination  4>  is  to  be  evaluated  at  |  =  $.  as  Indicated  by  the  subscript. 

On  the  other  hand,  the  cone  radius  rc  is  evaluated  at  R  »  R.  +  dR.  A  step 
by  step  application  of  Equations  (18)  and  (19)  yields  the  streamline  trace 
in  terms  of  its  coordinates  R  and  <)>.  The  pressure,  velocities,  etc.,  along 
the  streamline  can  then  be  determined  by  interpolation  with  respect  to  $ 
within  the  conical  solution  previously  obtained. 


The  final  ingredient  needed  for  the  boundary  layer  analysis  is  the  metric 
for  the  coordinate  normal  to  the  streamlines,  h2>  which  is  a  measure  of  the 

local  streamline  spreading.  This  can  be  determined  from  the  differential 
equat i on** 


d  An  h2 
— 


d  An  r 
_ c 

ds 


r^  V  3<ji 


(20) 


*fromS?h^UH!!riP^  j.here  distinguish  mesh  points  along  streamline  traces 
from  the  discrete  intervals  at  which  the  conical  solutions  are  obtained. 

actual  !y  an  approximate  version  of  Equation  (B-ll)  of  Refer- 
(  7  '  The  approximation  resides  in  the  fact  that  the  cone  radius  rc  ap- 
pearing  on  the  right  hand  side  should  actually  be  the  radius  of  the  equivalent 
body  of  revo  ution  as  defined  in  Appendix  B  of  Reference  (17).  For  the  re- 
has  !ha  °W  Jjone  h.a,f  angles  and  angles  of  attack  for  which  this  analysis 

*”c tnr  H?^e,°ped  kth,s  approximation  should  be  reasonable.  Note  that  the 

on  i  eT?nn  nWeen  r<:  3"d  the  fctual  curvaturc  can  be  estimated  to  be 

on  the  order  of  100  (I  -  cos  if  cos  6  ) . 

c 
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Again,  a  finite  difference  scheme  is  employed  to  integrate  Equation  (20) 
The  form  taken  is 


(h2)  =  exp  (Fj  +  F2) 

j  +  l 

Fj  =  In  [(rc)  (h2)  (rc)  '] 

j  +  l  j  j 

,e  i  3V.  .  3 V 

rr  .♦  ‘rW1.1 

C  J+l  C  r  J 

Initial  conditions  for  Equation  (20)  are  applied  at  R  *  dR  in  accord¬ 
ance  wi  th 

3 

h2  =  sin  6c  (dR) '  +  2  k 


(21a) 


(21b) 


(21c) 


where 

i  2_  3V/3P 
3  V  sin  8c 

Equation  (22)  gives  the  correct  limiting  behavior  for  zero  angle  of  attack 
(i  .a.^2  -►  r  )  and  is  consistent  with  the  known  solution  for  p  ■  0  as  de¬ 


veloped  by  Roshotko  (Reference  23) 


The  procedure  utilized  to  integrate  Equations  ( 1 8)  and  (21)  is  as 
fol lows : 


1.  Starting  at  a  generic  point  j  the  variable  R  is  incremented  by  an 

amount  dR  from  which  (r  )  =  (R.  +  dR)  sin  0  is  computed. 

j  +  l  J  C 

2.  Since  4* j  is  known,  the  Increment  dp  can  be  computed  from  Equation  (19) 


3.  The  increment  ds  then  follows  directly  from  Equation  (18) 
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4.  With  $  -  4> .  +  d$  interpolation  within  the  stored  conical  solution 

yields  the  value  of  (1/V  3V  /3$)  .  ,  needed  In  Equation  (21). 

4>  J  + 1 

5.  The  right  hand  side  of  Equation  (21a)  can  now  be  evaluated  yielding 
the  desired  value  of  (h,) 

J+l 

In  carrying  out  this  procedure  a  test  is  made  on  the  relative  change  of 
l/V  3V  /3d>  corresponding  to  a  choice  of  dR.  If  the  change  is  greater  than 
a  certain  tolerance  the  step  size  is  halved  and  the  integration  step  is  re¬ 
peated.  Thus,  the  step  size  for  integration  is  controlled  internally  by  the 
computer  logic  based  on  the  specified  tolerance.  Numerical  experiments  made 
using  the  computer  code  based  on  this  analysis  have  shown  that  a  relative 
change  of  5%  or  less  provides  sufficient  accuracy. 

The  computer  code  based  on  this  analysis  is  called  STRMIN  and  is  de¬ 
scribed  in  detail  in  the  user's  manual  which  is  available  from  the  AFFDL/FXG 
Numerical  results  obtained  using  this  code  and  an  assessment  of  its  ac¬ 
curacy  are  presented  in  Section  II  I A  celow. 

E.  Determination  of  Effective  Bo*/  Shape  -  The  approach  utilized  to  de¬ 
termine  the  effective  body  shape  involved  use  of  the  FFCS  described  in  Sec¬ 
tion  MB  above  in  conjunction  with  data  supplied  by  the  Air  Force,  Reference 
(24).  The  latter  consisted  of  surface  pressure  and  heat  transfer  measure¬ 
ments  around  the  periphery  of  a  6  half  angle  cone  approximately  36  inches 
long  and  pitot  and  total  temperature  distributions  throughout  the  shock  layer. 
These  measurements  were  obtained  at  a  free  stream  Mach  number  6  and  unit 

Reynolds  number  of  4.25  x  10^  per  inch  with  the  cone  at  angles  of  attack 
a  o  3 t  $ t  9  and  12  degrees.  These  tests  were  conducted  with  both  a  sharp 
cone  and  blunted  cone  (nose  radius  =  0.4  inches). 

The  effective  body  shapes  for  each  of  these  cases  was  established  by  an 
iterative  procedure  wherein  different  body  shapes  were  prescribed  as  input  to 
the  computer  code  until  a  satisfactory  match  of  the  calculated  surface  pressure 


. 


variation  with  the  measurements  was  obtained.  This  procedure  was  entirely 
successful  as  will  be  demonstrated  by  the  results  presented  in  the  next 
section.  Further  confirmation  of  the  adequacy  of  the  resulting  flow  field 
calculations  was  provided  by  comparison  of  the  calculated  and  measured 
pitot  pressure  distributions.  Again,  excellent  agreement  was  obtained. 

Results  of  these  comparisons  are  also  shown  in  the  next  section. 

Once  this  step  had  been  completed  it  was  possible  to  characterize  the 
effective  body  shape  dependence  on  angle  of  attack  by  relating  the  peripheral 
distribution  of  shear  layer  thickness  to  the  displacement  thickness  at  the 
primary  separation  point.  The  latter,  of  course,  was  estimated  by  means  of 
the  BLCS  described  in  Section  I  1C. 

At  this  point  in  the  development  it  was  necessary  to  generalize  the 

nominal  body  shapes  thus  established  to  account  for  the  effects  of  Mach 

* 

number  and  cone  half  angle  as  well  as  for  angles  of  attack  greater  than 
those  studied  above  (i.e.,  a/0  >  2).  This  generalization,  or  extrapolation 

would  be  established  by  utilization  of  other  data  which  is  available 
(cf:  Table  I)  as  well  as  the  pressure  correlations  which  were  developed 
under  the  current  ATL  effort.  These  pressure  correlations  are  an  essential 
ingredient  of  our  method  of  approach  and  are  discussed  in  the  next  section. 

F.  Leeward  Side  Pressure  Correlations  -  Since  the  number  of  geometric 
and  flight  parameters  which  define  the  flow  configuration  are  numerous 
(a,  0  ,  M^,  Re^)  and  their  range  of  interest  quite  extensive  it  was  recognized, 
at  the  outset  of  this  study,  that  detailed  analysis  of  data  for  all  conditions 
and  ranges  of  conditions  would  not  be  feasible.  Accordingly,  the  approach 
adopted  here  was  to  establish  the  basic  form  of  our  analytic  model  from  de¬ 
tailed  examination  of  a  limited  set  of  experimental  results  and  then  utilize 
the  remaining  data  to  "calibrate"  the  model  to  properly  account  for  variations 

in  M  ,  etc. 

00  ' 


■'Reynolds  number  effects  are  taken  into  account  automatically  via  the  dis 
placement  thickness  calculation  at  the  primary  separation  point. 
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1 nvest i gator 
and  Reference 

M 

00 

Re°°  X 
(millions) 

a (degrees) 

6c  (degrees) 

(s/r„) 

Anderson  (25) 

10.8 

.003  “  .04 

1-15 

10.0 

sharp 

Buck  (26) 

7.5,  ll.o,  14.0 

1.6  -  50.0 

5-20 

7.0 

25-80 

Feldhuhn  (3) 

5.07 

.5  -  5.0 

10-30 

5.0 

sharp 

Fitch  (27) 

10.2 

.86  -  8.5 

14.0 

7.1 

sharp 

Horstmann  (28) 

41.0* 

.04  -  .5 

5,  10,  15 

3.0 

sharp 

Rhudy  (24) 

6.0 

3.5  -  15.0 

3-12 

6.0 

sharp 

+ 

22-84 

Rainbird  (2) 

1.8,  4.25 

21.2,  43.4 

15.6,  22.8 

12.5 

sharp 

Stetson  (4) 

14.5 

.056  -  .80 

2-18 

5.6 

sharp 

+ 

8-92.5 

2.4-23.6 

Tracy  (1) 

7-9S 

.05  -  .4 

2-24 

10.0 

sharp 

Widhopf  (6) 

5.0 

1.0  -  50.0 

10-27.5 

9.0 

1 .05-5.25 

Yahalom  (5) 

2.72 

.35  -  .56 

5-30 

10,  15,  20 

sharp 

Zakkay  (29) 

6.0 

4.0  -  40.0 

7.5-22.5 

7.5 

sharp 

The  Air  Force  data  of  Reference  (24)  was  selected  as  the  primary  source 
for  the  model  development  since  it  provided  the  most  detailed  set  of  measure¬ 
ments  at  conditions  which  were  of  most  interest  from  a  practical  standpoint 
and  permitted  systematic  examination  of  the  effect  of  angle  of  attack.  The 
latter,  of  course,  is  the  essential  parameter  controlling  the  phenomenon. 

The  remaining  data  was  utilized  to  establish  a  correlation  which  could 
allow  estimates  to  be  made  of  the  pressure  ratios  p^/p^,  Pp/pw  and  Pg/PM  as 
defined  in  Figure  (4)  for  given  values  M^,  ot/6 c  and  Re^,  x,  the  latter  being 
a  Reynolds  number  based  on  free  stream  conditions  and  length  along  a  conical 
ray.  Although  some  of  the  data  sources  gave  anomalous  results,  a  satisfactory 
correlation  was  achieved  and  is  presented  in  the  next  section,  With  this  cor¬ 
relation  available  the  means  for  calibrating  the  analytic  model  for  any  com¬ 
bination  of  the  controlling  parameters  is  available. 

G.  Summary  -  ATL's  method  of  approach  may  be  summarized  as  follows: 

1.  The  separated  regions  of  high  shear  on  the  leeward  side  are 
represented  analytically  by  modifying  the  body  cross  sections 
in  an  appropriate  manner. 

2.  The  form  of  these  modifications  is  determined  by  an  iterative 
procedure  involving  the  Kutler  flow  field  program  and  selected 
pressure  data  -  a  match  of  calculated  and  measured  surface 
pressure  distributions  defines  the  desired  body  shapes. 

3.  Generalization  of  these  effective  body  shapes  for  practical 
ranges  of  interest  of  the  controlling  parameters  Is  established 
by  "calibration"  of  the  model  using  the  available  experimental 
data  in  the  form  of  pressure  correlations. 

4.  A  description  of  the  flow  field  throughout  the  shock  layer  Is 
obtained  by  means  of  the  Kutler  program  applied  to  a  body 
geometry  which  takes  on  the  effective  shapes  as  deduced  above. 
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SECTION  III 

NUMERICAL  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

A-  Streamline  Tracing  Scheme  -  Typical  results  obtained  with  the 
STRMLN  code  are  shown  in  Figure  (14).  These  results  correspond  to  one  of 
the  experimental  cases  of  Reference  (24)  so  that  the  requisite  pressure 
input  could  be  obtained  directly  from  measurements.  In  Figure  (15)  the 
pressure  distribution  determined  by  Equation  (8)  is  compared  with  these 
measurements.  As  can  be  seen  the  agreement  is  excellent  and  demonstrates 
the  applicability  of  the  form  selected  for  curve  fitting  of  these  pressure 
d i s tr i but i ons  . 

A  further  check  on  the  validity  of  this  formulation  Is  provided  by 
the  results  shown  in  Figure  (16).  Here  we  have  compared  the  flow  inclina¬ 
tion  i| ;  as  computed  by  the  STRMLN  code  with  that  predicted  by  the  FFCS. 

Further  indirect  evidence  of  the  satisfactory  performance  of  the  STRMLN 
is  provided  by  the  boundary  layer  results  presented  !n  the  next  section. 

B-  -B-our|dary  Layer  Computation  Scheme  -  To  verify  the  adequacy  o."  the 
selected  BLCS  a  series  of  calculations  were  carried  out  at  conditions  for 
which  experimental  data  were  available  for  comparison.  Representative  re¬ 
sults  are  shown  in  Figures  (17)  through  (21).  In  Figure  (17)  estimates  of 
Stanton  number  variation  along  the  several  streamlines  defined  in  Figure  (10) 
are  shown.  These  include  the  most  windward  streamline  (<j>  *>  0°)  as  well  as 

those  which  issue  (initially)  in  the  4>q  =  22*°,  47i°  and  60°  azimuthal  direc¬ 
tion;  i.e,:  at  zero  Incidence  these  streamlines  would  be  coincident  with 
the  cone  generators  defined  by  these  values  of  <f>. 

The  Stanton  number  estimates  are  compared  with  the  experimental  distribu¬ 
tions  deduced  from  the  Air  Force  data  (Reference  24)  in  Figure  (17)  and  indi¬ 
cate  excellent  agreement  in  both  the  laminar  and  fully  turbulent  regions.* 

’'Note  that  the  computer  code  descri  bed  in  Reference  (20)  has  both"  a  laminar  and 
turbulent  boundary  layer  capability  which  has  been  exercised  here.  This  pro¬ 
gram  also  has  a  built-in  transition  criterion  and  an  entropy  swallowing  option  as 
well.  For  the  present  purpose,  however,  the  program's  transition  option  was 
bypassed  and  transition  was  forced  to  occur  at  the  axial  station  implied  by  the 
data  as  indicated  in  these  figures.  These  options  have  been  included  in  the 
BLCS- 
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FIGURE  16.  COMPARISON  OF  PREDICTIONS  BY  STRMLN  AND  FFCS 
FOR  FLOW  INCLINATION 
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FIGURE  17  (CONTINUED) 
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FIGURE  18  STANTON  NUMBER  VS.  DISTANCE  ALONG  STREAMLINE  ON 
WINDWARD  PLANE  FOR  a  =  6° 


3- 


FIGURE  19  STATON  NUMBER  VS.  DISTANCE  ALONG  STREAMLINE  AT  WINDWARD 
PLANE  FOR  a  =  9° 
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Additional  comparisons  for  the  windward  streamline  at  angles  of  attack 
a  "  6°,  9°  are  presented  in  Figures  ( 1 8)  and  (19)  respectively.  Evidently, 
the  agreement  for  these  cases  is  as  accurate  as  for  the  a  *»  12°  results 
shown  in  Fi gure  (17)* 

A  demonstration  of  the  adequacy  of  the  BLCS  in  terms  of  Its  ability 
to  predict  integral  thickness  variation  around  yawed  cones  is  shown  in 
Figure  (20).  The  comparison  here  is  with  the  data  of  Rainbird  (Reference 
2).  To  obtain  our  theoretical  estimates  a  number  of  streamlines  were 
traced  at  the  appropriate  values  of  angle  of  attack  and  Mach  number  up  to 
an  axial  station  x  =  35  inches  corresponding  to  the  point  where  Rainbird's 
profile  data  and  i ntegral  th I cknesses  were  determined.  With  the  pressure 
distribution  and  spreading  parameter  h2  thus  determined,  the  BLCS  was 
exercised  yielding  the  variation  of  6  and  6  along  the  entire  streamline. 

The  desired  values  needed  for  comparison  with  the  experimental  results  were 
then  read  off  at  the  appropriate  axial  stations. 

The  results  are  shown  in  Figure  (20)  for  all  of  the  cases  examined  in 
Reference  (2).  As  can  be  seen,  the  agreement  for  both  integral  thicknesses 
is  satisfactory  up  to  peripheral  angles  of  $  %  120°. 

Additional  evidence  of  the  excellent  performance  of  the  BLCS  in  terms 
of  displacement  thickness  predictive  capability  is  demonstrated  in  Figure 
(21).  Here  we  have  compared  the  pressure  measurements  of  Reference  (2*0 
with  those  computed  by  means  of  the  FFCS  with  and  without  the  displacement 
thickness  effect  predicted  by  the  BLCS.  The  latter  distributions  are  shown 
in  Figure  (22).  It  is  interesting  to  note  that  despite  the  relatively  high 
Reynolds  numbers  which  are  involved  here  a  small  but  discernible  viscous 
interaction  effect  is  evident  even  on  the  windward  portions  of  the  cone. 

Finally,  in  Figures  (23)  and  (24),  results  are  presented  which  demon¬ 
strate,  indirectly,  the  satisfactory  performance  of  the  STRMLN  code  in  terms 
of  the  computation  of  the  streamline  spreading  parameter  hj.  Thus,  in  Figure 
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FIGURE  22  ESTIMATED  DISPLACEMENT  THICKNESS  DISTRIBUTION  OVER  SHARP  CONES 
AT  INCIDENCE 
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FIGURE  24  COMPARISON  OF  DISPLACEMENT  THICKNESS  AND  STANTON  NUMBER 
VARIATION  ALONG  THE  WINDWARD  PLANE  WITH  AND  WITHOUT  THE 
STREAMLINE  SPREADING  EFFECT. 

M  =  4.25;  9  =  12^°;  a  =  15.6° 


(23)  we  repeat  the  comparison  shown  earlier  in  Figure  (17a)  and  include  an 
additional  computation  due  to  the  BLCS  but  without  the  effect  of  streamline 
spreading.  That  is,  in  one  case  the  computation  was  carried  out  using  the 
h^  distribution  determined  by  the  STRMLN  while  in  the  second  case  the  actual 
geometric  body  radius  was  utilized  to  represent  the  h^  distribution.  As 
can  be  seen  in  Figure  (23)  the  Stanton  number  predictions  are  in  substantial 
disagreement  with  the  ~ieasu remen ts  when  the  body  radius  is  utilized,  par¬ 
ticularly  for  the  laminar  case,  while  the  prediction  based  on  the  h^  varia¬ 
tion  generated  by  the  STRMLN  provi des excel  lent  agreement  with  the  data. 

A  similar  result  is  shown  in  Figure  (24).  Here  we  have  compared  the 

variation  of  Stanton  number  and  displacement  thickness  on  the  windward  plane 

for  one  of  Rainbird's  cases  (M  =  4.25,  a  =  15-6°)  as  determined  with  and 

00 

without  the  streamline  "spreading"  effect.  It  is  apparent  that  the  effect 
of  streamline  spreading  on  the  Stanton  number  level  is  not  too  significant 
as  has  previously  been  pointed  out  by  Wionopf  (Reference  30).  Note,  for 
example,  that  at  x  *  35  the  percent  difference  between  the  two  values  of 
St  is  only  5%.  However,  the  corresponding  influence  on  displacement  thick¬ 
ness  is  far  from  negligible.  Indeed,  the  displacement  thickness  at  this 
station  would  be  over  predicted  by  about  60%  If  streamline  spreading  were 
neglected.  In  contrast,  with  the  spreading  effect  included  the  prediction 
is  within  approximately  10%  of  the  measured  value.  In  the  present  applica¬ 
tion,  of  course,  an  accurate  representation  of  integral  thicknesses  is  cru¬ 
cial.  Thus,  these  results  demonstrate  the  importance  of  utilizing  computa¬ 
tional  tools  which  properly  reflect  the  three  dimensional  character  of  both 
the  inviscid  and  viscous  flow  phenomena  which  is  of  interest. 

C.  Determination  of  Effective  Body  Shape  -  The  results  of  our  efforts 
to  deduce  appropriate  body  shapes  to  represent  the  viscous  layers  on  the  cones 
studied  in  Reference  (24)  are  summarized  in  Figures  (25)  through  (30).  Figure 
(25)  shows  the  final  selections  of  shear  layer  thickness  distributions.  The 
pressure  variation  computed  by  the  FFCS  for  these  effective  bodies  is  com¬ 
pared  with  the  experimental  data  in  Figure  (26).  Comparison  of  pitot  pro- 
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FIGURE  25  EFFECTIVE  BODY  SHAPES  UTILIZED  TO  COMPUTE  PRESSURE  DISTRIBUTIONS 
SHOWN  IN  FIGURE  (26) 


FIGURE  26  COMPARISON  OF  COMPUTED  SURFACE  PRESSURE  DISTRIBU 
TIONS  WITH  MEASUREMENTS  OF  REFERENCE  (24) 
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files  for  the  a  -  6°  and  12°  cases  is  shown  in  Figures  (27)  and  (28). 


For  the  a  «  12°  case  a  pitot  pressure  map  was  also  prepared  and  ts 
shown  in  Figure  (29a).  This  may  be  compared  with  the  corresponding  expert 
mental  map  shown  in  Figure  (29b)  as  provided  in  Reference  (24). 


Finally  In  Figure  (30)  there  is  presented  an  expanded  view  of  the  pitot 

maps  for  all  three  angles  of  attack  to  permit  a  detailed  comparison  of  the 

region  of  high  shear  (p.  /p„  <  0.7)  and  the  effective  body  thicknesses. 

t2  t2°° 

Recall  that  the  latter  were  determined  by  matching  the  peripheral  static 
pressure  distributions.  This  comparison  as  well  as  all  the  other  results  shown 
In  Figures  (25)  through  (29)  demonstrate  that  our  choices  of  effective  body 
shape  adequately  represent  the  viscous  phenomenon  occurring  on  the  cone  sur¬ 
face. 


D.  Pressure  Correlations  -  The  source  of  the  data  examined  for  this 

correlation  has  been  enumerated  in  Table  I.  Figure  (31)  summarizes  all  of  the 

results  obtained  from  these  sources  for  the  leeward  pressure  ratio  p^/p^  in 

terms  of  the  Reynolds  number  Re  .  No  effort  at  a  correlation  Is  intended 
'  °°,x 

here  but  the  data  is  presented  in  this  manner  to  provide  a  convenient  refer¬ 
ence  and  to  indicate  the  degree  of  scatter  which  exists. 

ATL's  effort  to  correlate  all  of  these  disparate  data  Involved  an  ex¬ 
amination  of  a  variety  of  parameters  with  the  initial  attempt  being  the  use 
of  the  correlation  suggested  by  Feldhuhn  (Reference  3).  These  and  other 
Intermediate  results  will  not  be  reported  here  but  can  be  found  detailed  in 
the  several  progress  reports  issued  under  the  current  contractual  effort. 

The  final  choice  of  correlation  schemes  involves  the  use  of  a  modified 
form  of  the  viscous  interaction  parameter  for  a  cone 

X  5  M  (C  /Re  )  Vs  i  n20 
C  “  »  °°,x  c 

The  modification  involves  a  combination  of  xc  with  the  relative  angle  of  in- 
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Ref.  (2) 


LEEWARD  PRESSURE  RATIO  VS  REYNOLDS  NUMBER -SUMMARY  OF  EXPERIMENTAL 


(25a) 


cidence  n  =  a/6^  according  to 

F4  -  n/xj 

Fk  =  (a/6  )  (M  /Re*  /sln20  )"*  (25b) 

'  c  <*>  «,x  c 

A  plot  of  all  of  the  data  presented  In  the  previous  figure  in  terms  of 
p  /p  vs  Fk  is  shown  in  Figure  (32).  Although  some  correlation  Is  apparent 
within  individual  sets  of  data  (see  especially  the  data  of  References  1  and 
26),  collapse  of  the  entire  spectrum  of  measurements  is  obviously  unsatis¬ 
factory.  On  the  other  hand,  if  the  correlation  scheme  is  restricted  to  the 
following  range  of  the  pertinent  variables  ** 


M  >5 

CO  — 

Re  >  5x10^ 

00  — 


a  satisfactory  correlation  is  obtained  as  shown  in  Figure  (33). 

The  restriction  with  respect  to  nose  bluntness  effect  is,  of  course, 
reasonable  and  simply  delineates  the  point  at  which  bluntness  effects  can  be 
neglected  and  the  flow  phenomenon  can  be  considered  self-similar  ( 1 . e . :  conical). 
The  restrictions  enumerated  above,  however,  do  serve  to  eliminate  several  data 
sets  obtained  on  sharp  cones.  Spec! fical  ly,  these  include  the  experimental  re¬ 
sults  of  Anderson  (Reference  25),  Yahalom  (Reference  5)  and  Rainbird  (Reference  2). 


’’'Note  that  the  selected  form  of  correlation  parameter  F4  effectively  sets  Coo 
equal  to  unity.  However,  the  effect  of  this  parameter  has  not  been  examined 
in  our  correlation  scheme. 

**Also  excluded  are  the  helium  results  of  Reference  (28). 
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>  5  x  10 


For  the  first  of  these,  no  rationale  for  their  exclusion  can  be  offered  at 
this  time  beyond  noting  the  very  low  Reynolds  number  involved  (order  of  10  ). 

The  rationale  for  rejecting  the  latter  two  sets  of  data  is  associated  with 
the  concept  of  the  "base  interference"  regime  postulated  by  Stetson  and  Ojdana 
(Reference  4).  It  is  believed  that  the  combination  of  low  Mach  numbers 
(1  •  8  <_  M^  <_  4.25)  and  thick  cone  angles  (10°  <_  6c  <_  20°)  for  which  the  in¬ 
vestigations  of  References  (2)  and  (5)  were  conducted  places  their  results  in 
the  aforementioned  category.  It  should  be  noted,  however,  that  the  restric¬ 
tions  associated  with  this  correlation  are  consistent  with  the  requirements 
and  objectives  of  the  subject  contract.  Accordingly,  the  selected  correlation 
shown  in  Figure  (33)  is  considered  to  be  satisfactory  for  the  present  purpose. 

To  completely  characterize  the  pressure  on  the  leeward  side  of  cones 
at  incidence,  a  way  of  estimating  the  pressure  ratios  Pp/p^,  end  ar'd 

their  corresponding  peripheral  locations  and  as  defined  in  f-lgure  (4) 
is  also  required.  For  evaluation  of  p^  and  an  analysis  similar  to  that 
described  above  for  p^  was  carried  out.  The  results  of  these  studies  led  U 
the  final  form  of  correlation  shown  in  Figure  (34).  According  to  this  cor¬ 
relation,  for  values  of  F4  %  0.7  the  minimum  pressure  occurs  at  the  leeward 
plane  of  symmetry;  i.e.:  the  pressure  decreases  monoton i ca 1 1 y  from  the  wind¬ 
ward  to  the  most  leeward  ray.  As  F4  continues  to  increase  a  minimum  :n  pres¬ 
sure  (p  )  occurs  outboard  from  the  p  =  180°  plane.  However,  for  values  of 

P  .  . 

F4  up  to  about  2,  the  pressure  subsequently  increases  monoton i ca 1 1 y  without 

an  inflection  point  to  the  leeward  value.  Finally,  for  F4  ^  2  a  second 
"minimum"  in  pressure  and/or  an  inflection  point  occurs  In  the  pressure  dis¬ 
tribution  (p$)  . 

The  available  data  for  the  peripheral  locations  p^  and  pg  are  presented 

in  Figure  (35)  plotted  against  the  relative  Incidence  n.  As  can  be  seen  these 

parameters  appear  to  be  essentially  independent  of  n  and  nose  bluntness  for  the 

ranges  considered.  On  tie  basis  of  these  results  it  is  concluded  that  suitable 

estimates  for  p  and  p  are 
P  5 
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FIGURE  35.  LOCUS  OF  FIRST  AND  SECOND  PRESSURE  MINIMUMS  AS  A  FUNCTION  OF 


Generalization  and  Calibration  of  the  Effective  Body  Shape  -  The 


selected  effective  body  shape  is  shown  schematically  in  Figure  (36).  Complete 
definition  of  this  geometry  requires  specification  of  a  total  of  seven  para- 

Vc  3 *«  V?  3'; 

meters;  viz:  6^,  6^,  <5^,  <S^  and  and  <j>s  where,  for  the  most  general  case, 


these  would  all  be  a  function  of  M  ,  a,  9  and  Re 

»  c  °°,x 


On  the  basis  of  our 


analysis  to  date  the  following  provisional  representation  is  employed.  Let 

n  =  a/0  .  Then 
c 


k_ (n) 

C  p 


ss  =  ks(n)  6p 


m  k.  (n)  6" 

L  L  p 


*C  *  F1(n) 


*$  "  p2(n) 


6n  =  G. (M  ,  n,  9  ,  Re  ) 
0  1  «  c  °°,x 


6  ■  G,(M  ,  n,  8  ,  Re  ) 

p  2  “  c  «,x 


Our  explicit  choices  for  the  functional  forms  appearing  in  Equation 
(27)  are  shown  in  Figure  (37).  The  discrete  symbols  which  appear  therein  re¬ 
present  the  actual  values  which  were  employed  In  the  FFCS  to  repeat  the  cal¬ 
culations  previously  made  with  the  body  shapes  shown  in  Figure  (25).  The  sur¬ 
face  and  pitot  pressure  distributions  which  resulted  were  und I st 1 ngui shabl e 
from  those  presented  in  Figures  (26),  (27)  and  (28). 


PROVISIONAL 

MODEL 


USED  IN 
CALCULATIONS 


FIGURE  37  VARIATION  OF  EFFECTIVE  BODY  COEFFICIENTS 
WITH  RELATIVE  ANGLE  OF  INCIDENCE 


The  parameters  4q  and  6^  represent  the  displacement  tnickrasses  at  the 
windward  plane  ($  ■  0)  and  at  the  o  *=  125  plane.  These  are  computed  case  by 
case  utilizing  the  BLCS  and  procedures  described  in  Section  !  s/B . 

We  have  indicated  above  that  this  represen ta t i on  is  provisional.  This 
follows  from  the  fact  that  the  coefficients  appearing  in  Equat ions (27)  have 
only  an  n  -  dependence  which  was  deduced  from  a  single  set  of  experimental 
results  obtained  at  a  fixed  value  of  M  and  0  and  for  a  relatively  limited 
range  of  the  parameter  n  (0.5  "  -  <  2.0).  At  this  time  the  additional  cal¬ 
culations  needed  to  establish  whether  or  not  the  representation  shown  in 
Figure  (37)  can  be  considered  universal  and  how  these  curves  can  be  extra¬ 
polated  to  larger  values  of  n  have  not  been  undertaken.  This  would  con¬ 
stitute  the  'calibrati on' '  phase  of  the  study.  Howe ve  r  ,  some  indication  of 
t.he  generality  of  our  provisional  model  has  been  obtained  and  is  demonstrated 
by  the  results  shown  in  Figure  (30).  Here  we  have  compared  the  maximum  thick- 

■jV 

nesses  (4^)  from  the  present  results  with  those  measured  by  Yahalom  (Reference 

5)  at  conditions  which  differ  markedly  from  those  of  Reference  (24).  The 

agreement  of  the  two  sets  of  data  is  remarkable  and  provides  some  justification 

for  considering  the  parameters  k  (n),  ..  etc.  independent  of  both  M  and  6  . 

c  «  c 
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SECTION  IV 
CONCLUDING  REMARKS 


The  methodology  described  In  the  previous  sections  is  capable  of  pro¬ 
viding  a  complete  description  of  the  flow  field  throughout  the  shock  layer 
surrounding  a  circular  cone  at  large  angle  of  attack.  The  method  Is  equally 
applicable  to  sharp  and  blunt  cones  provided  that,  in  the  latter  case,  the 
region  of  interest  lies  beyond  a  point  25  nose  radii  downstream  of  the 
stagnation  point. 

In  order  to  implement  the  method  the  combined  use  of  three  distinct 
analytic  tools  is  required.  Each  of  these  is  implemented  by  a  computer  code 
and  user's  manuals  for  each  of  these  codes  is  provided  in  Volume  II  of  this 
report . 

The  sequence  in  which  these  codes  are  utilized  is  as  follows: 

(a)  The  FFCS  is  applied  to  the  geometric  cone  at  the  selected 
angle  of  attack  to  provide  an  estimate  of  the  pressure  dis¬ 
tribution  around  the  cone  up  to  a  peripheral  angle  4>  ~  120°. 

This  information  is  needed  as  input  to  STRMLN. 

(b)  The  STRMLN  code  is  used  to  trace  a  series  of  streamlines  on 
the  surface  of  the  geometric  cone  up  to  b  -  120°. 

(c)  BLCS  is  used  to  estimate  the  boundary  layer  properties  at  b”0  and 
$-120°  based  on  the  results  obtained  from  step  (b)  above. 

(d)  The  boundary  layer  properties  obtained  in  step  (c) ,  when  input  to 
the  FFCS,  provide  the  means  for  estimating  the  effective  body 
shape  which  best  represents  the  viscous  effects  induced  by  the 
large  incidence.  Using  this  effective  geometry  the  FFCS  computes 
all  flow  properties,  including  the  surface  pressure  throughout 

the  shock  layer.  The  range  of  applicability  of  the  methodology  is 
restricted  by  limitations  inherent  in  the  FFCS  analysis  and  com- 
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puter  code.  Thus,  the  BLCS  and  STRMLN  analyses  and  codes  are 
considered  reliable  for  the  entire  range  of  conditions  which 
are  likely  to  be  of  practical  interest.  On  the  other  hand,  at 
the  present  time  the  FFCS  has  to  be  considered  provisional  since 
it  has  not  been  tested  over  a  sufficiently  wide  range  of  condi¬ 
tions.  This  is  particularly  true  of  the  general ized  effective 
bode  shapes  as  indicated  in  Section  I VE . 

In  AT L ' s  judgement,  the  methodology  can  be  expected  to  provide  reason¬ 
able  estimates  of  flow  field  characteristics  at  ^ach  numbers  on  the  order  of 
6  and  for  relative  angles  of  incidence  up  to  about  2.  For  conditions  which 
vary  to  any  substantial  degree  from  these  values  further  calibration  of  the 
method  would  be  required  before  it  can  be  considered  a  reliable  technique. 
ATL  believes  that  such  a  calibration  effort  is  feasible  and  that  it  would 
be  a  useful  undertaking  for  the  immediate  future. 
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